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Do double features in averaged pulsar profiles decipher the 
nature of their radio emission? 
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ABSTRACT 

An interesting paper has recently been published claiming that the long-sought Rosetta 
Stone needed to decipher the nature of pulsar radio emission has been finally identified 
as the double features in averaged pulsar profiles. The authors argue that highly sym- 
metric bifurcated features are produced by a split-fan beams of extraordinary-mode 
curvature radiation emitted by thin microscopic streams of magnetospheric plasma 
conducted by a very narrow bundle of magnetic field lines. We examined arguments 
leading to these intriguing conclusions and found a number of flaws. At least one of 
them is fatal, namely there is not enough available energy within such thin microscopic 
plasma streams. Using an elementary pulsar physics we show that if the stream is so 
thin that its emission can reveal the signatures of elementary radiation mechanism, 
then the energy deficit tends to be severe and reaches a few to several orders of mag- 
nitude (depending on the actual efficiency of converting the available kinetic energy 
of relativistic charged particles into the coherent radio emission). We are certain that 
the answer to the question contained in the title of this paper is definitely negative. 
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1 INTRODUCTION 

Almost half a century passed since the discovery of pul- 
sars, but yet no agreement has been reached concerning 
the actual mechanism for their radio emission (of genera- 
tion of their observed radio emission). The exceedingly high 
brightness temperature that can be deduced from the ob- 
served flux densities undoubtedly implies that the pulsar 
radiation must be emitted coherently. Generally, the pul- 
sar radio emission can be generated by means of either a 
maser-like or the coherent curvature mechanism (Ginzburg 
& Zheleznyakov 1975; Kazbegi, Machabeli & Melikidze 1991; 
Melikidze, Gil & Pataraya 2000, Paper I hereafter). Without 
any doubt this radiation is emitted in a strongly magnetized 
electron-positron plasma well inside the light cylinder. Once 
the waves are generated in the emission region, in the prop- 
agation region they naturally split into the ordinary and 
extraordinary waves corresponding to the normal modes of 
the strongly magnetized plasma (see e.g. Arons & Barnard 
1986; Lominadze et al. 1986). The ordinary waves are polar- 
ized in the plane of the wave vector k and the local magnetic 
field B and their electric field has a component along both 
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fe and B. Therefore, they interact strongly with plasma par- 
ticles (causing charge-separation along field lines), and thus 
encounter strong difficulty in escaping from the magneto- 
sphere. On the other hand, the extraordinary waves are lin- 
early polarized perpendicularly to the plane containing both 
fe and B vectors and thus they cannot separate charges 
along B. As a result the extraordinary mode can propa- 
gate freely through the magnetospheric plasma and escape 
to the interstellar medium. Many observational constrains 
on the emission altitude unambiguously suggest that the 
emitted radiation detaches from the magnetosheric plasma 
at altitudes r being less than 10% of the light cylinder ra- 
dius Rlc = 21x1 Pc (e.g. Blaskiewicz, Cordes & Wasserman 
1991; Kijak & Gil 1997). It is worth emphasizing that this 
also holds for the millisecond pulsars (see Figure 3 in Kijak 
& Gil 1998) 

Therefore, from a theoretical point of view the bulk of 
the observed pulsar radiation originates when the extraor- 
dinary plasma waves escape from the magnetosphere. There 
exists also strong observational indication that the extraor- 
dinary mode is dominant in pulsar radiation. Indeed, Lai, 
Chernoff & Cordes (2001) found that the Vela pulsar emits 
radio waves polarized predominantly in the direction per- 
pendicular to the plane of dipolar magnetic field lines. In 
fact, they were able to demonstrate convincingly that in 
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the fiducial phasqjthe radiation is polarized perpendicularly 
to the plane of the dipolar magnetic field lines. This argu- 
ment was extended by Gil, Lyubarsky & Melikidze (2004, 
Paper II hereafter ) to every phase within the pulse win- 
dow, using the fact that the mean position angle swing in 
this pulsar follows the rotating vector model (RVM here- 
after; Radhakrishnan & Cooke, 1969; Johnston, van Straten, 
Kramer & Bailes, 2001). This means that the radiation ob- 
served at a given longitude is polarized perpendicularly to 
the plane of dipolar magnetic field lines, along which the 
sources of this emission are moving. It is worth noting that 
the above conclusion concerns the average polarization. The 
analogous problem related to the instant emission observed 
in single pulses was recently discussed by Mitra, Gil & Me- 
likidze (2009, Paper III hereafter). 

The open question is which of the two possible mech- 
anisms of coherent radio emission: maser-like or curvature 
radiation, is responsible for generation of the observed pulsar 
radiation? Let us keep in mind that it must be the extraordi- 
nary mode (to escape freely from the magnetosphere) polar- 
ized perpendicularly to the planes of dipolar magnetic field 
lines (to satisfy observational polarization constrains). The 
most natural candidate is the coherent curvature radiation, 
as it is the only mechanism that distinguishes planes of mag- 
netic field lines (source trajectories). Recently Mitra, Gil & 
Melikidze (2009) found strong arguments to help distinguish 
between the alternative mechanisms. They analyzed highly 
polarized (nearly 100%) single pulses in a number of strong 
pulsars and argued that they can be produced only by the 
extraordinary mode excited in the magnetospheric plasma 
by means of the coherent curvature radiation. In fact, they 
found that position angle variations in subpulses precisely 
follow the RVM-like mean position angle swing. This is ex- 
actly what is expected to be produced by curvature radiation 
in a plasma, whose polarization is perpendicular to the mag- 
netic field line planes. The maser-like emission generates fast 
swings of instant position angle in subpulses, incompatible 
with the RVM (see Paper III for details). 

It would certainly be desirable to find additional and 
independent observational evidences to support the coher- 
ent curvature radiation being the mechanism for generation 
of pulsar radio emission. Recently, Dyks, Rudak & Rankin 
(2007; DRR07 hereafter) and Dyks, Rudak & Demorest 
(2010; DRDIO hereafter) claimed that the double symmet- 
rical features (called the bifurcated components; BFC here- 
after) and notches (absorption features) observed in aver- 
aged profiles of some pulsars (e.g. their Figs. 1 and 2) carry 
the crucial information able to decipher the nature of the 
observed radio emission. The idea was to fit these features 
with the elementary emission pattern of selected radiation 
mechanisms: the parallel acceleration beam and the curva- 
ture radiation beam. Although we do not believe that the 
actual pulsar radiation mechanism can be identified from the 
analysis of the properties of average waveforms alone, we do 
not intent to argue with the approach itself. Rather, we in- 
tend to verify the arguments and claims of DRDIO using 
different and independent methods. 



^ The fiducial phase corresponds to the fiducial plane, which 
contains both the rotation and the magnetic axes as well as the 
line-of-sight (that is the wave-vector k) 



Early attempts were unsuccessful and it turned out that 
the fit of the parallel acceleration beam model proposed by 
DRR07 "was the wrong idea". Most recently DRDIO ulti- 
mately gave up on the parallel acceleration mechanism and 
concentrated on the curvature radiation. They claimed that 
the observed pulsar radio emission was the coherent curva- 
ture radiation (more precisely the component polarized per- 
pendicularly to the planes of magnetic field lines). At first we 
were glad to see published conclusions strongly supporting 
our results obtained by means of completely independent ar- 
guments and methods (Papers I, II and III). Later we found 
a number of flaws in considerations of DRDIO. We decided 
to present and discuss these flaws in this paper, as we were 
afraid that misleading and incorrect arguments of DRDIO 
can rather harm the idea of the coherent curvature radia- 
tion as the mechanism of pulsar emission than promote it. 
Some readers may get an impression that the long standing 
problem of the physical mechanism of the pulsar emission 
has just been solved. Unfortunately this is not true and still 
more work is to be done in this field. 

First of all, the considerations of DRDIO are based on 
the single particle curvature radiation mechanism in vac- 
uum, while it is well known that pulsar radiation must 
be emitted coherently in the magnetized electron-positron 
plasma. Apparently, DRDIO assumed implicitly that a hy- 
pothetical coherence mechanism (which they did not spec- 
ify) would just reconstruct the single particle vacuum case. 
However this is not true. As we demonstrate in this paper 
the single particle vacuum model cannot even be used as 
a first approximation of coherent curvature radiation in a 
magnetized plasma. 



2 DOUBLE FEATURES AND NOTCHES 

As already mentioned, DRDIO ignored realities of the ex- 
citation of radio waves and their propagation in the pulsar 
magnetospheric plasma. All their arguments were based on 
the formalism of the single particle curvature radiation ex- 
cited and propagated in vacuum. They strongly concluded 
that the observed pulsar radiation was the coherent curva- 
ture radiation, although without any justiflcation for the co- 
herency part. They apparently believed that their arguments 
could be directly applied to the realistic pulsar environment. 
For example they state: "The bunching- induced coherency 
(e.g. due to the two-stream instability, Ruderman & Suther- 
land 1975; RS75 hereafter) seems to have less problems with 
the quasi- isotropic amplification of the non-coherent beam." 
This statement apparently mixes two independent problems, 
namely: creation of bunches and generation of coherent radi- 
ation. The bunches formed naturally by linear electrostatic 
waves (which DRDIO seem to be referring to) cannot pro- 
vide any emission, because the characteristic time-scale of 
such bunching is too short compared with the time-scale of 
curvature radiation (see Paper I for details). On the other 
hand, the RS75 mechanism corresponds to vacuum, where 
there would be no problem with generation of electromag- 
netic waves with both polarizations (provided they could 
be emitted). Thus, the problem DRDIO seems to be facing 
is how to dump the parallel mode of curvature radiation, 
which in vacuum is 7 times stronger than the perpendicular 
one (e.g. Jackson 1975). This happens naturally in magne- 
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tized plasma but at the same time properties of the escaping 
curvature radiation change their characteristics with respect 
to the vacuum case (Paper II). Thus, the latter cannot be 
used as a credible diagnostic tool to unravel the pulsar ra- 
diation mechanism. 

To the best of our knowledge the only efficient and phys- 
ically self-consistent mechanism of bunch creation in the 
pulsar magnetospheric plasma is the spark associated soli- 
ton model developed in Paper I. This model uses the two- 
stream instabilitjo exclusively to generate the longitudinal 
(non-electromagnetic) plasma waves. The actual bunching 
is caused by the nonlinear evolution of these plasma wave- 
packets and formation of charged relativistic solitons capa- 
ble of emitting coherent curvature radiation. The influence 
of the ambient magnetized plasma on both generation and 
propagation of this radiation was studied in Paper II, where 
it was shown that only the extraordinary mode (with polar- 
ization perpendicular to the planes of magnetic field lines) 
can escape from the pulsar magnetosphere. Following Paper 
II we will use the point-like approximation model of soliton 
bunches, making our arguments quite general, i.e. indepen- 
dent of the actual bunching mechanism. We will demon- 
strate that characteristic properties of the coherent curva- 
ture radiation in the pulsar magnetospheric plasma are quite 
different from that of the textbook vacuum case explored by 
DRDIO. We will examine an influence of the plasma on the 
frequency dependence of the bifurcation angle of the BFC. 
Additionally, we will check if there is enough kinematic en- 
ergy (the uppermost limit to the radiation energy) to power 
the BFC within the thin plasma stream model considered 
by DRDIO. 



2.1 Fitting the bifurcated component 

The most important feature considered by DRDIO is a 
bifurcated component (BFC) in the mean profile of PSR 
J1012-I-5307 (see their Fig.2). This feature has a high 
(although not perfect) degree of a mirror symmetry and 
DRDIO argue that it reflects a pure morphology of an ex- 
traordinary mode of curvature radiation Q Indeed, waves 
polarized perpendicularly to the plane of the charge mo- 
tion are not emitted along the instant velocity vector (see 
lower panel of Figure 3 in Paper II, where it is clearly seen 
that there is no curvature radiation emitted at and near 
the local magnetic field direction). DRDIO in their Section 
3.2.2 adopt a strong assumption concerning the origin of 
BFC feature, namely that "the emitter has a form of thin 
and elongated plasma stream that emits the curvature ra- 
diation mainly in the extraordinary (orthogonal) mode". 
However, their Eqs.(3) and (4) corresponding to the vacuum 
case contain two electromagnetic components, both of which 



^ It is worth emphasizing that this is the only plasma instabihty 
that can occur in the near pulsar magnetosphere (e.g. Asseo & 
Melikidze (1998). 

^ Strictly speaking they mean the component of curvature radi- 
ation polarized perpendicularly to the plane of source motion (i.e. 
in the plane of curved magnetic field line) , as neither ordinary nor 
extraordinary mode exists in vacuum. Both parallel and perpen- 
dicular polarization components of curvature radiation in vacuum 
have the feature of the extraordinary mode, which is defined by 
the absence of the electric field component along the wave- vector. 



should reach the observer. Moreover, in contrast to their as- 
sumption of missing the parallel component, its power is 
about 7 times higher than the power of the perpendicular 
mode. In order to dump the parallel component one needs 
a plasma environment, where the corresponding equations 
for the power of polarized radiation are significantly differ- 
ent from those in vacuum (see Paper II). DRDIO used the 
vacuum formalism to fit the curvature radiation to the BFC 
of PSR J1012+5307, while the proper approach would be to 
fit the coherent curvature radiation in plasma using a gen- 
eral formalism developed in Paper II (see next Section for 
details). Moreover, DRDIO fit both the parallel acceleration 
beam and curvature radiation beam models to the BFC pro- 
file and conclude that the latter is slightly better, although 
the former cannot be excluded. However, they found that 
only the curvature radiation model can match the data si- 
multaneously at two frequencies. In the next section we ex- 
amine this problem under the proper treatment including 
the ambient plasma infiuence. 



2.2 Frequency dependence of the bifurcation 
angle 

As mentioned above, DRDIO argue that the important prop- 
erty of BFC supporting the model of curvature radiation is 
the frequency dependence of the bifurcation angle (angu- 
lar separation between component peaks; see their Figs. 8 
and 9) Abfc oc i^~^g'^^ ■ Interestingly, the value of this expo- 
nent is close to 1/3 = 0.333, which indeed follows from the 
properties of the single particle curvature radiation in vac- 
uum (provided that Abfc = I/7 and r^obs oc 7^, e.g. Jackson 
1975). This would be an impressive feature supporting the 
curvature radiation model if one could be certain that the 
ambient plasma does not change the value of the bifurcation 
exponent. However, as we will demonstrate below, this is not 
the case. Moreover, even DRDIO admit that "the exponent 
of 1/3 is not ubiquitous among pulsars". In fact, as one can 
see in their Fig. 6 PSR J1012+5307 is the only pulsar having 
the bifurcation exponent close to 1/3 (it is also worth noting 
that the measured value 0.35 is not equal to 1/3 even within 
the error bars marked in the figure). In realistic theory of 
the coherent curvature radiation emitted and propagating 
in a pulsar plasma this exponent can differ from 1/3 and/or 
0.35. As an example we explored the formalism developed in 
Paper II, which corresponds to the general case of curvature 
radiation emitted by the point-like (smaller than the emitted 
wave-length) charged bunch/soliton moving relativistically 
along curved magnetic field lines. Only the extraordinary 
mode polarized perpendicularly to the planes of field lines 
can reach the observer. We calculated an opening angle be- 
tween the local magnetic field and the direction at which the 
maximum power of this mode is emitted in a plasma as a 
function of frequency (for illustration see the upper panel of 
Figure 3 in Paper II, where it is clearly seen that the extra- 
ordinary mode is missing towards the direction of the local 
magnetic field). The results are presented in Fig[T] where 
on the horizontal axis is the normalized observational fre- 
quency V jvcr (where v^i ~ 7.2 x 10^7'^/ p [Hz], and 7 is the 
Lorentz factor of the source of curvature radiation moving 
relativistically along the trajectory with the radius of cur- 
vature p), while on the vertical axis is the normalized angle 
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Figure 1. Frequency dependence of the opening angle of tlie 
extraordinary mode of colierent curvature radiation in the pulsar 
magnetospheric plasma (thick line). The formal power- low fit is 
represented by the dashed line and the 95% upper and lower limits 
are represented by dashed lines. 



ip/ifi^^ between the direction of maximum power emission 
and the local magnetic field vector (where ip^^ = I/7). 

Unlike the vacuum case, the bifurcation frequency de- 
pendence is not exactly a power-law like for the curvature 
radiation in a plasma (thick solid line in Fig. [T]). However, 
when the formal power-law fit was applied we obtained 
(p/(P;.j oc (://i/cr)~", where a = 0.45 ± 0.01. This value of 
bifurcation exponent is far from 1/3 that can be derived for 
the vacuum case. Moreover, it does not depend on the actual 
bunching mechanism (e.g. soliton model), that is any point- 
like emitter of curvature radiation will give the BFC expo- 
nent equal to 0.45 in the pulsar plasma. Also, the elementary 
pattern of the curvature emission in a plasma is different 
from that of emitted in vacuum. Moreover, the spectrum 
of the soliton curvature radiation differs from that of sin- 
gle particle even if both are calculated in vacuum (see Fig. 
4 in Paper I), not to mention the infiuence of the plasma 
environment (Paper II). The quality of fit of the curvature 
radiation pattern to the BFC profile in a plasma will cer- 
tainly be much worse than in the vacuum case. Thus, in our 
opinion there is no reason to believe that neither the formal 
fit itself nor the exponent value (0.35) in the frequency de- 
pendence of BFC found by DRDIO in PSR J1012-H5307 is 
the signature of the actual radiation mechanism, the curva- 
ture radiation in particular. If indeed there was a possibility 
to detect and resolve the elementary feature of the curvature 
radiation (from a single source or a number of sources fiow- 
ing along a narrow bundle of field lines and emitting in the 
same plane) then the bifurcation exponent should be about 
0.45 and not about 0.35. Moreover, as demonstrated in the 
next section there is not enough power within such a nar- 
row bundle to provide the observed BFC luminosity, which 
is even more fatal for the model proposed by DRDIO. It is 
worth emphasizing here that we do believe coherent curva- 
ture radiation to be the actual pulsar radiation mechanism, 
but for completely different reasons than those presented by 
DRDIO (see Papers I, II and III). We think that the BFC fea- 



ture in PSR J1012+5307 represents normal components in 
the complex mean profile of this pulsar and the measured bi- 
fircation/separation index is close to 1/3 by accident. There 
are many factors influencing a value of this index besides the 
intrinsic radiation mechanism, such as: radius-to-frequency 
mapping, radius of curvature varying across the emission 
region, size of the emission region, organization of the emis- 
sion region (conal versus patchy) and geometry of pulsar 
emission, (e.g. Gil & Krawczyk 1996, Gil et al. 2002). 



2.3 Energy considerations 

As already mentioned DRDIO used the textbook formalism 
of single particle curvature radiation in vacuum and implic- 
itly assumed that all properties of the observed coherent 
radiation will be exactly the same when considered prop- 
erly in the pulsar magnetospheric plasma (although they 
never expressed this explicitly). We conjecture that this as- 
sumption concerned the energy problem as well, which how- 
ever DRDIO did not consider at all. We can make simple 
estimates of the most upper limits of a possible emission 
power, and compare them with the observed radio luminosi- 
ties. The mean flux density from PSR J1012+5307 at 1.4 
GHz is about 3 mjy (ATNF database), which for the dis- 
tance d = 0.52 kpc translates into the radio luminosity Lr ~ 
6 X 10^^ erg s"^ (Eq.(3.14) in Lorimer & Kramer, 2005). 
Judging from the mean profiles shown in Figure 2 of DRDIO 
it is reasonable to assume that the BFC feature contains not 
less than about 10% of the total energy associated with the 
whole pulse profile. The same must hold for the emitted 
power, so we can assume that Lbfc ~ O.lLr = 6 x lO^'' 
erg s~^. The highest available energy source is determined 
by the spin down power Lsd ~ 4 x lO'^^P-is/P'^ erg s~^, 
which for PSR J1012+5307 with P = 0.0053 s and P_i5 = 
1.7 X X 10~^ gives Lsd ~ 4.7 x 10^^ erg s~^. We can now cal- 
culate the so-called pulsar kinematic luminosity Lkin, which 
is the power carried by charged particles accelerated within 
the pulsar inner gap (see Section 6.2 in Paper II). This lumi- 
nosity can be expressed as Lkin = J-pj-meC^najSpc erg s~^, 
where 7 < 5 x 10^ is the "primary" Lorentz factor of elec- 
trons (or positrons) leaving the acceleration region ("polar 
gap"), ncj = 1.4 X 10"(P_i5/P)°-^ cm"^ is the Goldreich- 
Jufian (1969) number density and Spc = 6.6 x 10**P"^ cm^ 
is the canonical polar cap surface area. For the parameters 
of this pulsar ncj = 8 x 10^ cm~^ and Spc = 1-3 x 10^^ cm^ 
and thus Ljtin < 10'^'^ erg s~^, which is few percent of L sd 
(as should be expected in general). The radiation efficiency 
77 of the observed radio emission referred to the kinematic 
pulsar luminosity is 77 = Lr/Lkin=6 x lO^'^/lO^^ = 6 x 10^^, 
which is quite typical for radio pulsars. Below we argue 
that such low efficiency should be also expected in conver- 
sion of the particle's kinematic luminosity into the power of 
coherent curvature radiation. Indeed, we can estimate the 
kinetic energy fiux of the coherently emitting bunches as 
Lkin ~ ■ymeC'^nGjS. If we assume that the entire Lkin is 
converted into the radio emission (i.e. Lr = Lkin), then the 
efficiency would be rj^j^ = Lkin /Lkin ~ l/lp^- For the typi- 
cal values of Lorentz factors 7 = 400 and 7 = 10® of the 
secondary and the primary plasma, respectively (see Papers 
I and II for details), we get r)Qi^ = 4 x 10"'* as a maximum 
possible efficiency estimation. More realistically, if only a 
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Figure 2. a. The sketch of the flux tube associated with the 
bifurcated component (BFC). b. The top view of the polar cap 
and the base of the BFC flux tube with a divergence Ax- 



part of Lkin is converted into Lr, then the efficiency is much 
lower, close to typical value measured in radio pulsars, i.e. 
'7cR ~ »? = 6 X 10-^ 

For a convenience of further considerations we will now 
introduce the surface density of kinematic luminosity, which 
can be defined as Lq — Lkin/5'pc. For the parameters of 
this pulsar Lq < 10^^ erg s^^ cm^^. According to DRDIO 
the BFC feature is emitted via curvature radiation of sources 
fiowing within a very narrow flux tube of magnetic field linos 
(see Figure 12 in their paper), although its actual cross- 
section is not specified. On one hand, this flux tube should 
be broad enough to carry much more than Lbfc = 6 x 10^^ 
erg/s in the particle flux. On the other hand, to reveal sig- 
natures of elementary emitters it should be narrow enough 
so that its divergence Ax = e/7 is much smaller than the 
opening angle of curvature radiation I/7. Thus, the auxil- 
iary dimensionless parameter £ << 1. Now we can roughly 
estimate the required cross-section of the flux tube satisfy- 
ing the above condition. We will use a projection along the 
dipolar field lines from the radiation region onto the stellar 
surface (see Fig. 2 for a schematic sketch). Near the edge 
of the polar cap we can write for the dimension of the base 
of the BFC flux tube and for its cross-section A = Ax'"pc 



and S — -KTp^Ax ~ 5'pcAx , respectively. Here rpc and Spc 
are the radius and surface area of the canonical polar cap, 
respectively. Now we can write that 



,2 2 

S^Spc^ = 1.3x 10" — , 



(1) 



and the associated kinematic luminosity carried along the 



a maximum emissivity at the frequency i/m = 2 x W^Y^/R 
Hz, where R > W^r^''' cm is the radius of curvature of 
the magnetic fleld lines and r is the emission altitude (see 
e.g. RS75). This frequency has to be about 1 GHz and 
thus, assuming reasonably that r > 3 x 10® cm, we obtain 
7 > 225. Therefore the kinematic luminosity of BFC flux 
tube L < 2.4 x lO^'^e^ erg/s. We still need to estimate the 
value of e which should be much less than unity. Let us note 
that even if we adopt as the most upper limit e ~ 0.3, we get 
the kinematic luminosity L < Lbfc = 3.6 x 10^® erg s^^. 
It is worth remembering that the extraordinary (perpendic- 
ular polarization) mode carries only 1/7-th part of the total 
emitted power, so the maximum efficiency of this mode is 
about 15%. This leads to the most upper limit L ~ 5 x 10^^ 
erg/s, much less than Lbfc = 6 x 10^® erg s~^. In prac- 
tice the actual efficiency is much smaller, but even assuming 
this absolutely unrealistic maximum efficiency there is not 
enough luminosity to power the BFC feature. 

The above estimate corresponds to unrealistic parame- 
ter values; 7 = 5 x 10®, 7 = 220 and e = 0.3. For more 
realistic values 7 = 10®, 7 = 400 and e = 0.1 (e.g. Paper 
I) one obtains L < lO^'* erg s^^ << Lbfc- Taking into ac- 
count that this kinematic power still does not include the 
efficiency of curvature emission we can see that the power 
deficit reaches several order of magnitudes. Once this effi- 
ciency is taken into account then one obtains 



Lr^< 



-^L = 2.5x lO^--' — 



BFC flux tube is L = LuS = 10 



^^£^7 



erg/s. In order to es- 



timate the value of 7 let us note that curvature radiation has 



(2) 

7pr 7 

The left-hand side of this equation represents the radio lu- 
minosity of the BFC feature, which should account for about 
Lbfc = 6 x 10^® erg s~^. For the first and the second set of 
values of e, 7 and 7 used above one obtains L^^ equal to 
0.03 and 6.2 times 10^° erg s~^, respectively, which is much 
less than Lbfc in both cases. Moreover, the inequality sign 
in the above equation is related to simplifying assumptions 
that we used, namely that charge density inside the coher- 
ently emitting bunch is equal to the local Goldreich- Julian 
(1969) value and the entire kinetic energy is converted into 
the radio emission, while in fact this is an upper limit. 

In summary, the model of the BFC emission presented 
by DRDIO has an energy deficit amounting to several or- 
ders of magnitude and it cannot be balanced by any means. 
This conclusion is independent of the actual pulsar radia- 
tion mechanism, i.e. it holds for any bunching mechanism 
(including the soliton model) in the pulsar plasma, the par- 
allel acceleration beam, maser-like emission, etc. The BFC 
feature can be emitted from the bundle of field lines base of 
which covers at least 10% of the polar cap area, once again 
irrespective of the actual radiation mechanism. 

We believe that BFC feature in PSR J1012-I-5307 is 
just a normal part of the multi-component average profile 
of this pulsar. What seems to distinguish it from the rest 
of the profile is its high degree of symmetry. However, such 
kind of symmetry is nothing extraordinary among the pulsar 
mean profiles. Moreover, one cannot exclude that the over- 
whelming symmetry of the BFC feature in PSR J1012-I-5307 
(see Figs. 3 and 4 in DRDIO) will disappear at different fre- 
quencies, which often happens in double component profile 
pulsars (some tendency to such change is visible in Figs. 8 
and 9 of DRDIO). Perhaps the best example of extremely 
symmetric profile is the case of PSR J063H-1036 (see fig. 2 
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in Weltevrede et al. 2008). Both the overall waveform and 
separation between the inner components in this profile is 
about the same as in the BFC of PSR J1012+5307. However, 
the excellent polarimetry available for the former pulsar in- 
dicates that its profile is emitted close to the fiducial plane^ 
and the separation/birfucation of several degrees of longi- 
tude is impossible for Lorentz factors 7 being of the order 
of 100 (necessary for the characteristic frequency of curva- 
ture radiation to be in the radio band). In PSR J1012-I-5307 
DRDIO solved this problem by postulating "the sightline 
cuts through plasma streams", but in PSR J0631-I-1036 the 
plasma streams are apparently cut more centrally. Other- 
wise, the double/birfurcated components look alike in both 
pulsars. 



3 CONCLUSIONS 

In the last paragraph of their paper DRDIO state "We con- 
clude that the long-sought Rosetta Stone needed to decipher 
the nature of pulsar radio emission has finally been iden- 
tified as double features in averaged pulse profiles". This 
statement, based on incomplete set of evidence, is simply 
not right. The BFC feature that DRDIO proclaimed the 
"Rosetta Stone" cannot be emitted in the way to play a role 
of the latter. There is not enough available energy source to 
power the features that could reveal signatures of the ele- 
mentary pulsar emission, no matter whether in single pulses 
or in average profiles. 

The huge energy deficit is the most serious problem of 
DRDIO model. It appears that the kinematic power begins 
to be balanced if the base of the bundle of magnetic field 
lines carrying the plasma stream associated with the BFC 
feature covers at least 10 % of the polar cap surface area, 
irrespective of the actual radiation mechanism. Of course, 
such stream is too wide to reveal the physical properties 
of the elementary emission, as there must be many sources 
of the coherent pulsar radiation distributed over its cross- 
section. 

Another major problem of DRDIO is the frequency de- 
pendence of bifurcation of the two components of BFC fea- 
ture (see their Figs. 8 and 9). DRDIO claim that this feature 
is produced by a split-fan beam of extraordinary-mode cur- 
vature radiation emitted along the sufficiently thin plasma 
stream. With no emissivity of the extraordinary-mode in 
the plane of field lines carrying this stream, the observed 
radiation should be bifurcated, with a bifurcation angle 
(angular distance between the two apparent components) 
being dependent on frequency. Using a very high quality 
observational data DRDIO found that Abfc oc v^^^, where 
a = 0.35 ± 0.01. They claimed that this was fully consistent 
with a low-frequency curvature radiation for which a — 1/3. 
First of all, the observed value is not equal to 1/3 even within 
error bars (see their Fig. 6), which is of course not a big 
problem. More serious problem is related to the fact that 
DRDIO used a textbook expression for the single particle 
curvature radiation emitted and propagated in vacuum. Al- 
though DRDIO never say it explicitly they implicitly sug- 
gest that the single particle vacuum curvature radiation is 
a very good approximation of the actual pulsar radiation 
mechanism (or in other words the influence of plasma on 
the generation and propagation of pulsar radiation can be 



neglected). While the single-particle radiation is not a bad 
model for a coherent radiation by a small point-like bunches, 
a vacuum approximation is absolutely not acceptable. We 
found that for the curvature radiation emitted and propa- 
gated in pulsar plasma (only the extraordinary-mode can 
leave the magnetosphere) the value of bifurcation exponent 
a — 0.45 ± 0.01. This should be the measured value of the 
bifurcation exponent if indeed it was possible to detect and 
resolve the elementary feature of the curvature radiation in 
the pulsar radio emission. Interestingly, there is one point in 
Fig.6 of DRDIO with a = 0.42±0.025 (J0437), but we think 
that this is by pure accident. Indeed, a "normal" separation 
exponents for components in the complex pulsar profiles can 
have any value between about 0.25 and 0.7, depending on 
different geometrical conditions (e.g. Gil & Krawczyk 1996, 
Gil et al. 2002 and Table 6 in DRDIO). We believe that the 
BFC feature is not different from the rest of the profile of 
PSR J1012-I-5307. Apparently, this is an almost aligned ro- 
tator and the polar cap is "seen" for almost an entire pulsar 
period. It would be extremely interesting to obtain a single 
pulse data, but this will probably be possible only with a 
future generation radiotelescopes. One important prediction 
is the following: if indeed the BFC feature in the profile of 
J1012-I-5307 represents the signatures of extraordinary mode 
of the curvature radiation then all single pulses beneath this 
component should also be bifurcated and look alike average 
emission. If it turns out not to be the case, then BFC feature 
is nothing extraordinary, just a normal macroscopic double 
component in multicomponent profile of this pulsar. 

The two major problems discussed above are fatal for 
the DRDIO model, especially the missing energy problem. 
Besides them there are few minor problems and we would 
like to mention two of them here. The natural consequence of 
DRDIO model is an elongated fan structure of pulsar beams. 
Such a structure is difficult to reconcile with the observed 
rates of occurrence of interpulses (IP), including double pole 
(DP) and single pole (SP) cases. Observational data show 
that in the entire pulsar sample there are about 2 and 1 per- 
cent of the former and the latter, respectively. Every statis- 
tical study of pulse profile widths that includes interpulses 
indicates that pulsar beams must be nearly circular (Gil 
& Han 1996, Kolonko, Gil & Maciesiak 2004, Weltevrede 
& Johnston 2008, Keith et al. 2010). Indeed, any signifi- 
cant deformation of beam circularity spoils the expected IP 
statistics. In particular, the elongated fan beams produce 
too many interpulses as compared with observations (e.g. 
Gil & Han 1996). 

Another minor problem is the following. The intrinsic 
beam- width of the BFC feature is less than one degree, while 
the observational bifurcation is about 8 degrees of longitude. 
To achieve the apparent broadening by a factor of 10 or 
so DRDIO must use extreme values of geometrical observa- 
tional angles (small cut angles) and/or substellar radii of 
curvature of magnetic field lines. This doesn't seem very 
likely, although such geometrical situation cannot be ex- 
cluded. 

We did not touched the phenomenon of notches at all, 
which seem to be bifurcated absorption features occurring 
in a few pulsars. According to DRDIO the physics of notches 
is identical to that of BFC and the only difference is geomet- 
rical in nature. Namely, the double absorption features are 
produced by beam of the extraordinary mode of curvature 
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radiation, when it is eclipsed by the thin plasma stream. We 
understand that this beam is analogous to the one associ- 
ated with the BFC feature discussed by DRDIO and thus 
the proposed origin of notches is subject to the same criti- 
cism as that of BFC related emission. However, the BFC is 
not visible and the only role the thin plasma stream plays is 
to eclipse a normal pulsar radio emission. This can happen 
due to geometrical (the thin beam emission misses the line- 
of-sight) or energetic (the beam is not producing the coher- 
ent radio emission) reasons. With this latter interpretation 
the origin of notches proposed by DRDIO can be perhaps 
considered as a viable model. However, it cannot be done 
without problems. One problem is how to pin firmly a thin 
plasma stream to the polar cap surface. The other, more se- 
rious problem is that the extraordinary mode of curvature 
radiation is not supposed to be absorbed by any kind of 
magnetospheric plasma. 

The criticism of DRDIO model presented in this pa- 
per does not change our opinion that the pulsar observed 
radiation is really the extraordinary mode of the coherent 
curvature radiation emitted and propagated in the magne- 
tospheric plasma (Papers I, H and HI). One is tempted to 
say that DRDIO came to the right conclusion for the wrong 
reasons. 
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